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Abstract 

We show that, given an arbitrary shift, the lapse N can be chosen so that the 
extrinsic curvature K of the space shces with metric g in arbitrary coordinates of 
a solution of Einstein's equations satisfies a quasi-linear wave equation. We give a 
geometric first order symmetric hyperbolic system verified in vacuum by g, K and A^. 
We show that one can also obtain a quasi-linear wave equation for K by requiring N 
to satisfy at each time an elliptic equation which fixes the value of the mean extrinsic 
curvature of the space slices. 
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SYSTEMES GEOMETRIQUES BIEN POSES POUR LES EQUATIONS D'EINSTEIN. 



Resume Nous donnons deux conditions differentes sur le lapse pour que la courbure ex- 
trinseque K des sections d'espace d'une solution des equations d' Einstein satisfasse une 
equation d'onde quasi-lineaire, le shift etant arbitraire. 

Version frangaise abregee 

Nous considerons les equations d'Einstein sur un espace temps de dimension 4, de topologie 
M X R dont nous ecrivons la metrique 

dans le corepere donne par (x* coordonnees locales sur M, t E R) 

00 = dt, 9' = dx' + (3'dt. 

Nous definissons I'operateur sur les tenseurs d'espace dependant du temps par {C est la 
derivee de Lie, Oq la derivee pfaffienne) 

do = d/dt — Cp, alors que do = d/dt — P^d/dx\ 

La courbure extrinseque K des sections d'espace M x {t} est telle que 

doQij = -2NKij. 

Nous utilisons la decomposition 3+1 du tenseur de Ricci de I'espace temps pour obtenir 
I'identite, (2.1) donnee au §2, generalisation de celle donnee dans [2] dans le cas /3 = 0, ou 
on a pose H — Kl et {ij) — ij + ji. Cette identite donne une equation d'ondes pour K dans 
les cas suivants: 

1. est choisi tel que dgN + N'^H = 0. 

2. N est choisi tel que H — h, une fonction donnee. 

2 



Nous donnons dans le cas 1 un systeme symetrique hyperbolique du premier ordre satisfait 
par les espaces temps einsteiniens. 

Introduction 

We examine the Cauchy problemi for general relativity as the time history of the geometry 
of a spacelike hypersurface. Constraints on the initial data, 'g, the space metric, and K, 
the extrinsic curvature, can be posed and solved as an elliptic system by known methods. 
However, the equations of motion for g and K, which are essentially the Arnowitt-Deser- 
Misner canonical equations, do not, despite their usefulness, manifest mathematically the 
physical propagation of gravitational effects along the light cone. These equations, of course, 
contain gauge effects and cannot, therefore, yield directly a physical wave equation, that is, 
a hyperbolic system with suitable characteristics. (For a review of the relevant geometry, 
see, for example, El.) 

In this paper we give two different methods for obtaining exact nonlinear physical wave 
equations from the equations of motion of g and K. One of these, which relies on a "har- 
monic" time-slicing condition, gives equations of motion completely equivalent to a first- 
order symmetric hyperbolic system with only physically appropriate characteristics. We 
construct this system explicitly. Among the propagated quantities is, in effect, the Riemann 
curvature. The space coordinates and shift vector are arbitrary; and, in this sense, the 
system is gauge-invariant. 

Our gauge invariant nonlinear hyperbolic system, being exact and always on the physical 
light cone, is well suited for use in a number of problems that now confront gravity theorists. 
These include large-scale computations of astrophysically significant processes (such as black 
hole collisions) that require efficient stable numerical integration, extraction of gravitational 
radiation with arbitrarily high accuracy from Cauchy data, gauge-invariant perturbation 
and approximation methods, and posing boundary conditions compatibly with the causal 
structure of spacetime. 



3 



1. NOTATIONS AND 3+1 DECOMPOSITION 



We consider a manifold V of dimension 4, which has the topology M x R, with M a C°°, 
3 dimensional manifold. We denote hy {x,t) G M x R a point of V. When we take 
local coordinates they will always be adapted to the product structure: x\ i — 1,2,3 are 
coordinates on M and x^ = t & R. We consider on y a pseudo-riemannian metric g of 
lorentzian signature ( — h ++). It induces a properly riemannian metric cjt on each space 
submanifold Mt = M x {t}. We denote by /? and N respectively the shift and the lapse of 
the foliation. We define a coframe by 



dt. 



dx' + (5'dt. 



The corresponding Pfaff derivatives da are 
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In this coframe the metric reads 



ds" = gaf.e'^e^ = -N\ey + g.^e'e^. 



We define the operator do on t dependent space tensors by {Cp is the Lie derivative with 
respect to (3) 



The extrinsic curvature K of the space manifold is 




The Ricci curvature of space time admits the 3+1 decomposition, with H — K^h 



R 



. . ^ -N-'doK,j + HKij - 2Ki„,K"'j - N-'WjdiN + R, 




N-\VhK^j-djH), 



R% = -{N-^V%N - KijK'^ + N-%H), 



where an overbar denotes a quantity relative to the space metric g = {gij). 
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2. SECOND ORDER EQUATION FOR K 

In the formula above Rij, like the right hand side giving its decomposition, is a t dependent 
space tensor, the projection on space of the Ricci tensor of g. We compute its do derivative. 
First we compute d^Rij. The infinitesimal variation of the Ricci curvature corresponding to 
an infinitesimal 5g variation of the space metric is, with {ij) = ij + ji (no factor of 1/2) 

This expression applies to {d/dt)Rij with 6gij = {d/dt)gij and to CpRij with 5gij = C^gij. 
Therefore, using the relation between dogij and Kij, we obtain 

-2NR\,^K^h - NR^(,K,{^. 
We now use the expressions for R^i and Rij to obtain the identity 
a, = doR^, + V(,(iV2i?^.)0) = 

-do{N~^VjdiN) - NVidjH 
-V^,{Kj)hd^N) - 2NR^,,^K^h - NR^f^.K,)^ + HV,d,N. 

3. HYPERBOLIC SYSTEM FOR g, K, N 

We shall eliminate at the same time the third derivatives of and the second derivatives 
of H as follows (cf. this elimination in the case of zero shift ini). We compute 

We find at once 

d - - dN 1 - d - d 

-VAN ^ V,a.— - -/'(V(,-^?,), - Vi-g.,)dkN, 
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and we find an analogous formula when the operator d/dt is replaced by Cp. Finally 

d^V,d,N = VjOAn + {V(,(iVir,)0 - Vi{NK,,)}d'N. 

The third order terms in N and second order terms in H in the second order equation for 
K can therefore be written in terms of 

dj = N-^Vjdi{doN + N'^H). 

We satisfy the condition Cij = by requiring N to satisfy the differential equation (note 
that doN = doN) 

(N') doN + N^H = 0. 

Using the expression for H, the equation (A^') reads 

dolog{N{detg)-'/'} = 0. 

We find, generalizing the result obtained ini in the case of zero shift, that the general solution 
of this equation is 

N = a-^det{gy/'^, Boa = 0. 

The second equation is a linear differential equation for the scalar density a, depending only 
on p. The above choice of is called an algebraic gauge. A possible choice if /? does not 
depend on t is to take a independent of t and such that = 0. In practice, one may 
simply regard (N') as a differential equation to be solved simultaneously with (K') below. 

Taking into account the equation (N') satsified by A^ we see that the Einstein equations 
Rai3 = Pa(3 imply the wave equation for K 

(K') NUK.j = + e,,-, 

where we now set 
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+{2Wi^,Kj)i)d^N - 2HN~^diNdjN - 2d(^iNdj)H - 3dhNV^Kij 

An immediate consequence is the following theorem. 

Theorem In algebraic gauge and with the definition Kij = —{2N)~^dogij, the system 
Qij — Qij = is a quasi- diagonal third order system for gij with principal operator Odo- 
This system is hyperbolic if g is properly riemannian and N"^ > 0. (In practice, one can 
regard the definition of Kij as an equation (g') for d^gij, and then consider the differential 
equations ft^'), (K') and (g') to be solved simultaneously.) 

The usual local in time existence theorem for a solution of the Cauchy problem results 
from the Leray theory of hyperbolic systemsB. 

4. VERIFICATION OF THE ORIGINAL EINSTEIN EQUATIONS 

We set 

Suppose we have solved the equations 

a, = d,R,j + V(,(Ar2i?0)) = (4.1) 
that is, equivalently, doing a few computations, 

VoS., + V(.(Ar2E°)) - g,,Vh{N^TP'') + = 0, (4.2) 

where fij is linear and homogeneous in T^ap- 

We suppose that the stress energy tensor T^p = p^p — ^gapP satisfies the conservation 
laws VaT'^^ = 0. We deduce then from the Bianchi identities the equations 
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V'^E,;^ = 0. (4.3) 

Prom (4.2) and (4.3) we obtain for E^^j a quasi diagonal third order system with principal 
operator the hyperbolic operator doO. The vanishing for i = of E^/j results from the 
constraints satisfied by the original initial data 

E°|t=o = 0, E°°|t=o = 0, 

and the determination of 9oi^|t=o by the equation 

Ejj|t=o — 0. 

We then deduce from (4.2) and (4.3) the vanishing for t = of the derivatives of order < 2 
of Ea^. It results therefore, from the uniqueness theorem of Leray[4] for hyperbolic systems, 
that a solution of (4.1) with initial data satisfying the constraints satisfies the full Einstein 
equations E^^^ = 0. 

5. FIRST ORDER SYSTEM (vacuum) 

The preceding results can be extended without major change to dimensions greater than 
4. We will show that in dimension 4 a solution of the vacuum Einstein equations, together 
with the gauge condition {N'), satisfies a first order symmetric system, hyperbolic if g is 
properly Riemannian and > 0. Such a system could be useful to establish a priori 
estimates relevant to global problems, as well as to the physical applications mentioned in 
the Introduction. 

We have obtained for the unknowns g, K, and N the equations 

dogij = -2NKij, (5.1) 
do{N^) = -2NH{N''), (5.2) 
NOKij = Qij. (5.3) 
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To obtain a first order system we take as additional unknowns: 

N~'^doKij = Lij, VhKij = Mhij, di log N = ai 
N'^dodi log = aoi, aij = VidjlogN. 



We take as equation (|0 



The equation ( |5.3| ) gives 



doL,^ - NV^Mh^j = N{HU, - Q,,). (5.4) 

In three space dimensions, the Riemann tensor is a linear function of the Ricci tensor. Using 
the equation Rij = to express Rij, we write Qij as a polynomial in the unknowns and in 
g^^ . {Lij is the essential piece of the spacetime Riemann tensor R°ioj-) 

We have the following lemma proved for instance by using the commutativity d^di = dido 
on components of tensors 

Lemma For an arbitrary covariant vector Ui we have 

and an analogous formula for tensors with additional terms for each index. 
Using this lemma we see that Lij and M^j must satisfy the equation 



doMk,j - NWkU, = Ni^akL,, + M^^Kj^i + Ki^iMj^k' - Kl^^M' j^k (5.5) 
+Ki(^i{K^j)ak + aj)K\ - a^Kj)k) 



On the other hand, (|5.2| ) implies 

doai = -N{Hai + Mik^) (5.6) 
while (|5.2| ) and the lemma yield 

doaji - NVjaio = Nai (M(.,j)^ - M\j + a^^Kj-^ - a^Kij) + Najaio. (5.7) 



Now we deduce the value of doaoi from (|5.2| ) and the Einstein equation 
R% = -{N-'^V''VhN - KijK'^ + N-^doH} = 0. 
Indeed these two equations imply 

Hence, by differentiation and use of the Ricci formula and the definitions of ahi and a^i, we 
obtain 

doaoi - NV\k^ = Niy-R\au + ai{H^ - 2KkiK''' + 2a''ak + 2a\) (5.8) 
+2aka\ + HMik'' - 2K''M,m 



We use again the equation Rij = in ( p. 81 ) to replace Ry by its values in terms of the 
unknowns. 

The right hand sides of the equations Q, Q), (U), (U), (U), and 



( ^.8|) are polynomial in the unknowns and g^^; they do not depend on their derivatives. We 
eliminate the first derivatives of g appearing on the left hand sides for instance by introducing 
on M an a priori given metric e which may depend on t but is such that 

doCij = 0. 

We denote by D the covariant derivative in the metric e. We deduce then from ( |5.1|) the 
following equation for Ghij = DhQij- 

doGhi, = -2N{ahK,, + Mu, + ^'"m*^.)-}- (5-9) 
The tensor S*, the difference of the connections of g and e, is given by 

ij — 2^ ~ ^hij)- 

We thus obtain a first order differential system, equivalent to a 
symmetric system, hyperbolic if A^^ > and gij is properly riemannian. Its characteristics 
are the physical light cone and the "time" axis do- 
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6. MIXED HYPERBOLIC ELLIPTIC SYSTEM FOR K, g, N WHEN H IS GIVEN 



Another procedure to reduce the second order equation for K obtained above to a quasi 
diagonal system with principal part the wave operator is to replace in the term VidjH the 
mean curvature if by an a priori given function h] this method was used in [3], with h = 0, 
in the asymptotically euclidean case. With this replacement the equations = take 
the form 



where Pij depends on K and its first derivatives and on g, N and OqN together with their 
space derivatives of order < 2. When and Qij are known the above equations together 
with (5.1) are again a third order quasi-diagonal system for g, hyperbolic if A" > and g 
properly riemannian. On the other hand, the equation i?"o = P°o together with H — h 
implies the equation 



This equation is an elliptic equation for N when ^, K and p are known. Note that 
for energy sources satisfying the "strong" energy condition we have — p°o > as well as 
|Arp = KijK^^ > 0. These properties are important for the solution of the elliptic equation. 

The hyperbolic system that we have constructed (for given A^ and 6) has local in time 
solutions for Cauchy data in local Sobolev spaces, hke the one obtained in algebraic gauge. 
But now A^ is determined by an elliptic equation which has to be solved globally on each 
space slice at time t. An iteration procedure leading to a local in time solution has been 
used in the vacuum asymptotically flat case with h — [3] . We give below a theorem which 
applies to the case of a compact M and which can also be proven by iteration. (Further 
theorems with h ^ on asymptotically flat spaces can also be given.) 

Theorem. Let (M, e) be a given smooth compact riemannian manifold. Let there be 
given on M X 7, 7 = [0, T], a "pure space" smooth vector field /3 and a function h such that 



(6.1) 



= {KijK'^ - p%)N = -doh. 



(6.2) 
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{Hk is the usual Sobolev space relative to (M, e)) 

he f] C^-^{I,Hk), doh>0, doh^O, /i(0, .) > 0, or /i(0, .) < 0, h ^ 0. 

2<k<3 

There exists an interval J = [0,£], i < T such that the system (6.1), with Qij — 0, (5.1), 
(6.2) with p°o — 0, has one and only one solution on M x J, 

n C'-'UH,), Kef] C^-\J,Hk), Are n C'-\J,H,^,) 

l<k<3 0<k<2 0<k<2 

with N > and g uniformly equivalent to e, taking the initial data 

9ij{0, •) = lij e Hs, Kij{0, .) = kij e H2 

if 7 is a properly riemannian metric uniformly equivalent to e and k.k ^ 0. 

It can be proved by using the Bianchi identities that if the initial data satisfy the con- 
straints the constructed solution satisfies the original Einstein equations. 

J.W.Y. acknowledges support from the National Science Foundation of the U.S.A., grants 
PHY-9413207 and ASC/PHY 93-18152 (DARPA supplemented). 
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